Mesoscopic charge fluctuations in the Coulomb blockade regime 
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We study mesoscopic fluctuations of the differential capacitance of a dot connected to a lead 
by a single-mode quantum channel with adjustable conductance G. We show that the amplitude 
of the fluctuations reaches maximum at a partially opened channel, G < e 2 /nh. Parametric corre- 
lations of fluctuations at different values of the applied voltage and magnetic field can be studied 
experimentally, and we find the corresponding correlation functions for the differential capacitance. 

PACS numbers: 73.23.-b, 73.23.Hk, 73.40. Gk 
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Charge of a conducting dot weakly connected to leads 
is quantized in the units of electron charge e. When the 
potential induced on the dot is varied, its charge varies 
in a step-like manner. (The potential can be tuned by 
means of voltage V g applied to a gate, which is coupled to 
the dot only electrostatically.) The step-like behavior of 
the dot charge Q(V g ), i.e., Coulomb blockade, was clearly 
demonstrated in experiments utilizing metallic islands 
jjj. The charge quantization is progressively smeared 
with the increase of the conductance Go of a junction 
connecting dot with a lead. Theory Q predicts that de- 
viations of Q from the quantized values are proportional 
to conductance, \Q(V g ) — me\ oc Go at Go <C G q . When 
Go is of the order of the quantum unit G q = e 2 /ttTl, the 
function Q(Vg) depends on the detailed properties of the 
junction For a single-mode junction, the Coulomb 

blockade vanishes when the reflectionless propagation of 
an electron mode through the dot-lead channel is almost 
reached, with the periodic part of Q(V g ) approaching zero 
approximately as G q — Go, see ||. Single- mode junc- 
tion allowing for ballistic propagation of electrons is a 
reasonable model for a semiconductor quantum dot de- 
vice. In such a device, Go is tunable, unlike in the litho- 
graphically prepared metallic samples [Q. This allowed 
several groups to attack the challenging problem of quan- 
tum charge fluctuations in the Coulomb blockade regime. 
Suppression of the Coulomb blockade effect at Go — >• G q 
was demonstrated experimentally by Waugh et al. [0] and 
Molenkamp et al. @ on double-dot devices. In very re- 
cent experiments, Berman et al. measured the charge 
smearing in a single dot. 

Mesoscopic fluctuations of the charge Q and corre- 
spondingly of the differential capacitance Cdiff of a dot 
reflect the randomness of the electron states within the 
dot. Gopar et. al. jj| employed the Random Phase Ap- 
proximation to consider the charge dynamics of the dot 
connected to a reservoir by a reflectionless (Go = G q ) 
one mode channel. In this approach, the mesoscopic 
fluctuations are associated with the fluctuations of the 
density of states on the Fermi level. Later it was shown 
JlO| that even in this regime the Coulomb interaction 
leads to the formation of a collective state involving one- 



electron states in a wide energy strip around the Fermi 
level (width of the strip is of the order of the charging 
energy Ec). As the result, the mesoscopic fluctuations 
of the capacitance at Go = G q are larger than predicted 
in Ref. H , and the scale of the correlation magnetic field 
is controlled by the charging energy rather than by one 
electron level spacing. 

In this Letter, we will study the capacitance fluctua- 
tions in the whole region < Go < G q . We will show that 
the mesoscopic fluctuations are determined by the fluctu- 
ations of a large number (~ Ec/ A, where A is the one- 
electron level spacing) of single electron wave functions 
rather than by the state at the Fermi level. Concentrat- 
ing on the low temperature limit T < A, we show that 
the dependence of the variance (SC 2 ^) on conductance 
Go is non-monotonous, with the maximum of fluctuations 
reached when the channel is partially open, Gq < G q . We 
will find explicit analytic expressions for the correlation 
functions characterizing the evolution of the fluctuating 
differential capacitance with the magnetic field. 

We account for the Coulomb interaction by adding an 
extra term to the Hamiltonian: 



H = H( S +E c {h-Nf 



(1) 



Here Hq is the Hamiltonian of the dot-lead system with- 
out Coulomb interaction, and the term Ec{n — M) 2 ac- 
counts for the charging energy of the dot, h = J2 P <4> c p * s 
the total number of electrons in the dot. The gate volt- 
age V g is represented by the conventional dimensionless 
parameter Af = V g C g /e, with C g being the gate capac- 
itance. The differential capacitance of the dot Cdiff is 
related to the ground state energy Eq of the Hamilto- 
nian (@) by 
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Here Q is the dot charge, en, averaged over the ground 
state of the system, G is the capacitance of the isolated 
dot; we have introduced a lever-arm coefficient G jC g into 
the definition of Cdiff ■ 

The strength of tunneling through the dot-lead junc- 
tion is characterized by its dimensionless conductance 
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g = Go/G q . To demonstrate the nonmonotonous be- 
havior of the amplitude of the differential capacitance 
fluctuations, we consider analytically the limits of weak 
tunneling (j< 1) and of almost ballistic propagation of 
electrons in the junction (1 — g -C 1). 

Weak tunneling. The case of the weak tunneling, 
g <C 1, is adequately described in the tunneling Hamilto- 
nian formalism, 

k p 
+ J2(t kp a{c p + t* kp c f p a k ) , (3) 

k.p 

where £fc and £ p are the one-electron energy spectra of 
the lead and dot respectively. Tunneling matrix elements 
here are related to g by g = viv d (\tkp\ 2 } /2n 2 . Here- 
after (...) denotes the ensemble averaging, and vi and 
Vd are the averaged densities of states in the lead and 
dot respectively. In the extreme g — > 0, the charge of 
the dot in the ground state is integer, Q — me. The 
half-integer values of the dimensionless voltage M corre- 
spond to the degeneracy points of the charging energy in 
Eq. (Q). If V g is tuned away from a degeneracy point 
(to — 1/2 < TV < to + 1/2), it takes a finite energy 
2Ec{N — m ± 1/2) to increase or decrease the number 
of electrons in the dot by one with respect to the ground 
state. 

The perturbation theory in the tunneling Hamiltonian 
(|J) can be used to calculate the dot capacitance (||), if the 
system is not too close to one of the charge degeneracy 
points 



g? ln(l/U) < 1 



(4) 



Here U = min{\N — to ± 1/2|}. The lowest-order correc- 
tion to the ground state energy is 



E osc {N) 



fJ 
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dtg d (t;Af)Gi(-t) 



(5) 



where Qd and Gi are the one-electron Green functions in 
the dot and the lead, determined by the Hamiltonians 
Hd and Hi respectively, and R c is the point of the dot- 
lead contact. The density of states per unit area in the 
dot and in the lead are denoted here as Vdo and vio- The 
dependence of Gd on the gate voltage is given by 

Gd(t,N) = g d0 (t) [0( t ) e -*(m+l/2-AO*ot 



where Gdo{t) is t ne one-electron Green function in the 
dot determined by the Hamiltonian Hdo — J2 P ^p c P c P - 
The exponential factors account for the change in the 
Coulomb energy of the system when an extra electron or 
hole tunnels into the dot. 



It follows from Eqs. (g), (||), and (g) that the fluc- 
tuations of the differential capacitance SCdiff are caused 
by the fluctuations of the one-electron Green function 
SGdo = Gdo — (Gdo) m the dot at the point of the contact. 
Then the statistical properties of 5C d ig can be related to 
the correlation function (SGd(ti; Bi)5Gd(t2', -B 2 ))|r c . The 
applied magnetic field B affects the Green function Gdo 
thus affecting the capacitance fluctuations. 

For a conducting dot, the correlation function for SGdo 
is conventionally expressed in terms of the classical prob- 
abilities, diffuson V and Cooperon C JO ]: 

($Gdo{ui; B^SGdoi^; B 2 )}\ Rc = ^v d aQ{r^\^2) 

x [v^d^-^D + c^-^d^-^D] | Rc . (7) 

For times exceeding the time of electron propagation 
across the dot of area S or, equivalently, at energies much 
less than the Thouless energy for a closed dot, Et, the 
probability density of an electron is homogeneously dis- 
tributed over the dot, and 
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(8) 



Here Q± = E T ± S 2 )/$ ) , and $ = 2irhc/e is 

the flux quantum. Differentiating Eq. (^|) twice over AT 
and extracting the random part, we come to the expres- 
sion that relates fluctuations of the dot differential ca- 
pacitance SCdiff to the fluctuations of the one-electron 
Green function in the dot 5Gd- The correlation functions 
for SCdiff are thus related to those for SGd, which are 
given by Eqs. (|)-(|). 

First we consider the variation of the dot capacitance 
with the magnetic field at the same gate voltage. The 
corresponding correlation function is 



(SCdiff (JV, B^SCdiffW, B 2 )) 
_ C 2 g 2 A 
~ 2Air i E c W 



B C (U) 



(9) 



with B ± = \B 1 ±B 2 \, U = min{|A/"- to± 1/2|}, and 



6 f 00 x 2 1 



V 2 + 2y 



y + i 



21n(l + y) 



The function M(x) cx 1/x 2 at large values of the argu- 
ment x 3> 1, so that the correlation function decays with 
the increase of B± over the value B c . At zero magnetic 
field, the variance of the capacitance is 



(SC 2 



diffl 
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The characteristic scale for the variation of the corre- 
lation function Eq. (O) with the magnetic field is [12] 



b c = 5° UL 

S V Et 



(10) 
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The meaning of this formula for the correlation field B c 
can be clarified if we recall the relation of the charge fluc- 
tuations to the electron Green function. In the semiclas- 
sical approximation, Q{t) oc e lS ^\ where S is the action 
calculated along the trajectory. The dwelling time of an 
extra electron in the dot is t d ~ Ti/U . Over this time, 
the trajectory of a diffusing electron would cover an area 
Dtd ~ TiD/U, if there were no bounds imposed by the 
dot size. The ratio rid = hD/SU ~ Et/U determines 
the typical number of "winds" for the electron trajectory 
inside the dot. Because of the random direction of the 
winds, the total flux threading the electron trajectory is 
<!> = BSy/nd, and the correction to the action due to 
magnetic field is Sb — 27r$/$o- Equating Sb to 2n, we 
obtain the formula for B c . 

Relatively short excitation lifetime, td <C fi/A, leads 
to correlations of the capacitance fluctuations in different 
Coulomb blockade valleys. In a valley, ~ Ti/Atd levels of 
spatial quantization closest to the Fermi level in the dot 
contribute to the charge fluctuations. Each subsequent 
peak of the Coulomb blockade corresponds to one more 
level occupied. Therefore, the fluctuations of capacitance 
at different valleys are correlated, unless the distance n 
between them exceeds Ec / A, so that an entirely new set 
of discrete levels contributes to the fluctuations. Calcu- 
lation of the corresponding correlation function yields: 



(6C diff {U)6C diff {U + n)) 
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(11) 



It is to be mentioned that Eqs. (||) and (^) are valid 
only if condition (Q) is satisfied. If N is close enough to 
one of the charge degeneracy points, so that (Q) is vio- 
lated (but U is still much greater than A/Ec), the finite- 
order perturbation theory is insufficient for calculation of 
the dot capacitance. Utilizing the renormalization group 
technique , we found that the capacitance fluctuations 
in this case are described by the same Eq. (g) with the 
renormalized junction conductance 



9 = 9 



cos 



25, 1 



(12) 



With this modification, the condition for the applica- 
bility of Eq. (§) is relaxed from (kU) to a weaker one, 
t/2 \n(l/U) < 1. This new condition is imposed by the 
requirement that the renormalized conductance g must 
be still much less than unity. 

Strong tunneling. As it was pointed out in the intro- 
duction, the effects of the Coulomb blockade become less 
pronounced as the dimensionless conductance of the dot- 
lead junction, g, approaches unity. In the limit A = 0, 
Flensberg H and Matveev || studied the case of a dot 



connected to the lead by a constriction with small re- 
flection amplitude r = y/1 — g. Since the characteris- 
tic energy Ec is much smaller than the Fermi energy, 
one can linearize the spectrum of one-dimensional elec- 
trons in the constriction. The electrons in the constric- 
tion are thus divided into left- and right-moving species, 
^{x) = e- lkFX ^ L (x) + e lkFX ^ R {x). In [§, the full Hamil- 
tonian (|l|) describing the electrons in the dot and lead 
connected to each other, is transformed to an effective 
one-dimensional form and expressed in terms ipL, ipR- In 
the limit A — > and in the absence of backscattering in 
the constriction, there is no mixing between the L and R 
states: 

Hm = iv F / dx j ip L d x tp L - ip R d x ip R \ 

J — oo 

+ E C (^J oIx-.^l+^r--^ ; (13) 

here vf is the Fermi velocity, and : . . . : stands for the 
normal ordering. The two terms in the Hamiltonian ( |l3| ) 
account for the kinetic energy of the electrons and charg- 
ing energy, respectively. Both scattering from a barrier in 
the constriction, and scattering from within a finite-size 
dot, create terms oc (ip R ip L + h.c.) in the Hamiltonian, 
and in the corresponding action [[l0| . The backscattering 
terms lead to essentially non-perturbative corrections to 
the ground state energy Eq. In the work of Matveev ||, 
the problem with a single scatterer in the constriction 
described by 



H bs = v F 



(14) 



was solved, and Eq was found. According to 0, the regu- 
lar Coulomb blockade oscillations in the energy are given 
by E osc {Af) - 2£cM 2 cos(27rA/')ln [l/|r| 2 cos 2 (ttAA)] at 

Backscattering from inside the cavity (A ^ 0) gener- 
ates a random contribution to Eq, superimposed on the 
regular oscillations. Once the reflection in the channel is 
accounted for non-perturbatively, the random backscat- 
tering from the walls of the dot can be studied by means 
of the lowest-order perturbation theory |hJ : 



v 3 r°° 

5E dt 5g d0 (t)G LR (-t;Af)l= R 

no J-oo 



(15) 



where Glr is a component of the electron Green func- 
tion in a ID infinite wire, described by the Hamiltonian 
H = H 1D + H bs [see Eqs. @, ©], ' 

G LR (t) - (^ R m L (0) + 4(t)iM0)>l*=o ■ (16) 

Function Glr describes the dynamics of backscattering 
of an individual electron from the dot. An electron wave 
in the channel carries charge and spin degrees of freedom, 
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and there are two time scales involved in the backscat- 
tering process. The shorter one, t p ~ h/Ec, corresponds 
to the reflection of the charge by the Coulomb barrier 
(resulting from the charging energy). The longer one, t a , 
describes the reflection of the spin component, and at fi- 
nite r is of the order of l/E osc (Af). At the intermediate 
time scale, 1/Ec < \t\ < t a , the charge and spin com- 
ponents of the backscattered electron are separated, and 
the Green function manifests a distinctively non-Fermi- 
liquid behavior, 



G LR {t;Af) 



;ttJV / 7 \t\ 



irv F t 1/ IT 3 t a (J\f) 



In 



1*1 



(17) 



Here 7 = e c with C = 0.577 being the Euler con- 
stant. The electron backscattering process is completed 
at t ~ £ CT , and at larger time scales the system exhibits the 
standard Fermi- liquid behavior with GLFt(t) ~ 1/t. The 
result @ 

can be obtained in the first-order perturbation 
theory in H os , with the low-energy cut-off E ~ E osc (N). 

Equations ^ and ( |l5| ) enable us to relate capacitance 
fluctuations SCdiff to one-electron properties of the sys- 
tem described by the functions SQdo and Glr- The major 
contribution to the integral ( |l5| ) is given by the interval 
1/Ec < \t\ < t a , where approximation ( |i~7| ) is applicable. 
The correlation function of the capacitance fluctuations 
is then calculated with the help of Eqs. (]?]) and (|J): 

{8C m (M,B 1 )5C diff (M,B 2 )) = lgc 2 A| r |2 

1 \ , / 1 



x cos 2nAf cos 2irM In 



|2 pnc2 . 
2 



In 



r| 2 cos 2 nA4 



K = l-g 



(18) 



with r max = max{|rcos7rA/'|, \r C0S7rA^|} and B c = 
($0 / S) \J2Ec I Et ■ The variation of the correlation func- 
tion with the magnetic field is described by the function 
P{B 2 /B 2 c ,r 2 max ). In the domain B < B c , 



P [—-r 2 

J^2 5 max 



1 3 



In- 



In max < 1 ; 



B 2 



B 2 
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At magnetic fields larger than the correlation field B c , the 
correlations rapidly decay, P (B 2 / ' B 2 ;r 2 ma ^) cx (B c /B) 4 . 
The perturbation theory Eq. ( Jl5| ) is applicable if the 



channel conductance is not too close to Go, 



» 



A/Ec- At a larger conductance, the amplitude of fluc- 
tuations saturates at SC - C(A/E C ) ln 2 {E c /A) @. 

Tunneling of arbitrary strength. The behavior of the 
capacitance fluctuations in the whole domain < Go < 
G q (with G q = e 2 /irh) follows from Eqs. (g) and (ph. 



At Go -C G q the fluctuations increase with conduc- 
tance, SCdiff /C ~ {Go/G q )yJ A/Ec- For an almost 
open channel (G q — Go G q ) the characteristic am- 
plitude of the capacitance fluctuations is SCdiff /C ~ 



[1 - (G /Gq)]y/A/E c , up to a In 3 factor. Thus the 
mesoscopic capacitance fluctuations reach maximum at 
a partial transmission of the channel connecting dot and 
lead. The parameter y/A/Ec for typical semiconductor 
dots G3] is ~ 0.2, so that the relative magnitude of these 
fluctuations is appreciable. 

In conclusion, we found mesoscopic fluctuations of 
charge and differential capacitance of a partially opened 
quantum dot. The amplitude of fluctuations reaches 
maximum at the intermediate junction conductance, 
Go ^5 e 2 /irfr. We also found the correlation function 
of these fluctuations at different magnetic fields and gate 
voltages. In the experiments ||-[8) no special attention 
was paid to the mesoscopic fluctuations. However, the 
parameters of quantum dots in J6| || are in the right do- 
main for studying mesoscopic fluctuations of the charge 
and of the differential capacitance. This can be accom- 
plished by observing the Coulomb blockade with a mag- 
netic field applied, which causes random variations of the 
charge. 
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